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ABSTRACT 


To better understand the early galactic production of Li, Be, and B by cosmic ray 
spallation and fusion reactions, the dependence of these production rates on cosmic ray 
models and model parameters is examined. The sensitivity of elemental and isotopic 
production to the cosmic ray pathlength magnitude and energy dependence, source 
spectrum, spallation kinematics, and cross section uncertainties is studied. Changes in 
these model features, particularly those features related to confinement, are shown to 
alter the Be- and B-versus-Fe slopes from a naive quadratic relation. The implications 
of our results for the diffuse 7-ray background are examined, and the role of chemical 
evolution and its relation to our results is noted. It is also noted that the unmeasured 
high energy behavior of a + a fusion can lead to effects as large as a factor of 2 in the 
resultant yields. Future data should enable Population II Li, Be, and B abundances to 
constrain cosmic ray models for the early Galaxy. 
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1 Introduction 


The Population I abundances of 6 Li, Be, and B (LiBeB) have been thought for some 
time to have their origin in spallation and fusion processes between cosmic ray and 
interstellar medium (ISM) nuclei (see, e.g., Reeves, Folwer, k Hoyle 1970; Meneguzzi 
Audouze, k Reeves 1971; Walker, Mathews, k Viola 1985). In the past few years these’ 
elements have been sought extreme Population II dwarfs, stars well known to exhibit 
the Spite plateau” in lithium (Spite k Spite 1982), which is understood to indicate the 
primordial Li abundance (Walker, Steigman, Schramm, Olive, k Kang 1991). Recently 
these same Pop II stars have also been shown to additionally contain Be, (Rebolo et al. 
1988a; Ryan et al. 1990, 1992; Gilmore et al. 1992a, 1992b; Boesgaard k King 1993) B 
(Duncan Lambert, k Lemke 1992) and most recently 6 Li (Smith, Lambert, k Nissen 
1992). These abundances provide important clues about the early galaxy, and cosmic 
rays have been considered the most likely production mechanism for LiBeB in these 

r!! (Steigman & Walker 1992 ( SW ); Prantzos, Casse, k Vangioni-Flam 1992 
(PCV); Walker et al. 1993 (WSSOF); Steigman et al. 1993 (SFOSW)). 

If indeed the extreme population II LiBeB arise from cosmic ray interactions, the 
study of their isotopic abundances opens an important window on astrophysics. In 
principle, we may be able to gain insight on early cosmic rays, as well as early star 
format ion rates and chemical evolution (e.g. PCV; Fields, Schramm, k Truran 1993; 
bilk k Schramm 1992). In addition, associated with these early cosmic ray events is 
an appreciable gamma-ray flux which would contribute (perhaps significantly) to the 
present diffuse gamma-ray background (Silk k Schramm 1992; Prantzos k Casse 1993). 
The early cosmic ray scenario is also of importance to cosmology. One may use the results 
of a galactic cosmic-ray (GCR) spallation and fusion model to help infer primordial 
Li from the observed Li abundance (Olive k Schramm 1993). One may also use the 
GCR model’s cosmic ray spallation information on the Li isotopes, compared with their 
observed abundances, to deduce the amount of possible stellar depletion (SFOSW) 

Our previous work (WS, WSSOF, SFOSW) was an attempt at a relatively model- 
independent approach. Without assuming a specific model (of cosmic ray or galactic 
chemical evolution) we concentrated on testing the consistency of standard GCR models 
with the recent LiBeB observations and with primordial nucleosynthesis where possible 
by examining elemental and isotopic abundance ratios. In another approach, Gilmore 
et al. (1992b), Feltzing k Gustaffson (1993), Prantzos (1993), Prantzos k Casse (1993), 
Pagel (1993), and particularly PCV have chosen to adopt a particular detailed model of 
early galactic chemical and cosmic ray evolution to examine its predictions and compare 
with observations. 

In this paper we will examine in detail the dependence of the LiBeB abundances 
produced by GCR nucleosynthesis on the uncertainties of, or allowed variations in, the 
cosmic-ray model. As such, we will discuss the uncertainties in cosmic ray models. After 
a brief review of the data, we present some model options in section 3. Our assumptions 
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regarding the evolution of the abundances of a, C, N, and 0 are discussed in section 
4. We address the issue of the Be and B slopes versus [Fe/H] in section 5. The cosmic- 
ray spectrum and the confinement of cosmic rays will be discussed in section 6. We 
also explore implications of these models on 7-ray production in section 7. We draw 
conclusions in section 8. 


2 LiBeB Abundance Data 

To make this work as self-contained as possible, we show in table I LiBeB isotopic 
abundances observed in Pop II halo dwarf stars. We list those stars in which at least two 
light elements have been observed as we will for the most part be primarily interested 
in elemental or, even better, isotopic ratios. The notation we employ below and use 
throughout the paper is that [X/H] represents the log abundance relative to the solar 
abundance, namely log(X/H) - log(X/H) @ and [X] = 12 + logpC/H). In the table, the 
iron abundance represents an unweighted “world” average. For the other abundances, a 
weighted average is given. The 6 Li abundance was taken from Smith et al. (1993). The 
7 Li abundances were taken from Spite k Spite (1982,1986); Spite et al. (1984); Hobbs 
and Duncan (1987); Rebolo et al. (1988b); Hobbs k Thorburn (1991); and Pilachowski 
et al. (1993). The 9 Be abundances were taken from Rebolo et al. (1988a); Ryan et al. 
(1992); Gilmore et al. (1992a, b); Molaro, Castelli, k Pasquini (1993); and Boesgaard k 
King (1993). Finally, the boron abundances were taken from Duncan et al. (1992). 

The ratios of 6 Li to 7 Li, Li to Be and of B to Be are the observed ratios. Because, 
for Pop II, the dominant contribution to the 7 Li abundance comes from primordial 
nucleosynthesis rather than GCR nucleosynthesis, a certain degree of caution is necessary 
when comparing the first two of these ratios to the predictions we discuss below. On the 
other hand because there is no appreciable big bang source for either Be or B (Thomas 
et al., 1993), the ratio of these two may be compared (unless there is an additional 
primary source for n B as discussed in Dearborn et al. (1988), Woosley et al. (1990), and 
Olive et al. (1993)). 

Because of the importance of the B/Be ratio, we note that the values given in the 
table for Be and the ratio B/Be, are averages over observations by several groups. 
These observations themselves show some spread which may be significant. For the 
star HD 19445, beryllium upper limits were obtained by Rebolo et al. (1988a), giving 
[Be] < 0.3, and B/Be > 1.3, Ryan et al. (1990) found the upper limit [Be] < —0.3 and 
hence B/Be > 5 which should be compared with the value of 3.5 in the table which 
represents the only positive identification of Be in this star by Boesgaard and King 
(1993) and is slightly discrepant with the upper limit of Ryan et al. For HD140283, we 
have measurements of Be by three groups: [Be] = -1.25 ± 0.4 from Ryan et al. (1992) 
giving B/Be = 14 ± 14; [Be] = -0.97 ± 0.25 from Gilmore et al. (1992) giving B/Be 
= 7 ± 5; [Be] < -0.90 from Molaro, Castelli, k Pasquini (1993), giving B/Be > 6; and 
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P e L~ from Boes gaard and King (1993) giving B/Be = 5 ± 3. Finally 

for HD201891, [Be] = 0.4 ± 0.4 from Rebolo et al. (1988a) giving B/Be = 20 ± 26 and 
[Be] - 0.67 ± 0.1 from Boesgaard and King (1993) giving B/Be = 11 ± 10. As one can 
see there appears to be a wide range in values (and uncertainties) in this ratio, and the 
values given my be completely dominated by systematics which are poorly accounted 
tor in the stated error. In particular, ratios are most useful if constant surface tempera- 
ture and surface gravity assumptions are used in the specific elemental determinations. 
Unfortunately, this is not yet the case. 


Indeed it is the uncertainties in these measurements and in the averages of these 
measurements which are themselves uncertain. Namely, the treatment of systematic 
errors is far from systematic. The conversion of line strengthes to abundances involves 
inferences on the surface temperature and surface gravity of the star. Most observational 
determinations have been made using differents sets of inputs. Though one can ascribe 
some uncertainty to chosen values of these inputs, it is not always clear to what extent 
these systematic errors have been incorporated into the quoted so-called statistical error, 
and different authors make divergent assumptions on the uncertainty of their assumed 
stellar parameters. Furthermore any average will surely underestimate the true error 
because of the mistreatment of systematic errors. Thus it is our feeling that B /Be and 
Li/Be ratios are extremely uncertain. 


Table 1. Observed Pop II abundances of LiBeB isotopes 


Star [Fe/Hj 


HD16031 

-1.9 

HD19445 

-2.1 

HD84937 

-2.2 

HD94028 

-1.6 

HD132475 

-1.6 

HD134169 

-1.2 

HD140283 

-2.6 

HD160617 

-1.9 

HD189558 

-1.3 

HD194598 

-1.4 

HD201891 

-1.3 

BD23°3912 

-1.5 


Li [Be] 

2.03 ±0.2 -0.37 ±0.25 

2.07 ±0.07 -0.14 ±0.1 

2.11 ±0.07 -0.85 ±0.19 
2.10 ±0.08 0.44 ±0.1 

2.05 ±0.09 0.60 ±0.3 

2.20 ±0.08 0.71 ±0.11 

2.08 ± 0.06 -0.87 ±0.11 

2.22 ±0.12 -0.47 ±0.18 

2.01 ±0.12 0.89 ±0.33 

2.00 ±0.2 0.37 ±0.12 

1.97 ±0.07 0.65 ±0.1 

2.37 ±0.08 0.30 ±0.4 


[B] c Li/ Y Li 

0.4 ±0.2 

0.05 ±0.02 

- 0.1 ± 0.2 
1.7 ±0.4 


Li/Be* B/Be* 
251 ± 185 

162 ±46 3.5 ±1.8 

912 ±425 
46 ±13 
28 ±20 
31 ±10 

891 ±262 6 ± 3 

490 ± 244 
13± 11 
43 ±23 

21 ± 6 11 ± 11 
117 ± 110 


* Ratios are extremely uncertain due to inadequate treatment of systematics. 
Formal errors on ratios are underestimates. 
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Cosmic Ray Models for LiBeB Spallation Pro- 
duction 


Currently studied models for LiBeB production are based on the work of Reeves, Fowler 
and Hoyle (1970) and the more detailed follow-up work of Meneguzzi, Audouze, and 
Keeves (1971). They describe cosmic ray propagation via the simple leaky box model, 
lhis assumes a spatially homogeneous distribution of sources, cosmic rays, and inter- 
stellar material. Within this model, the propagation equation is 


Ma 7 , d(b A N A ) 
dt A + dT 



( 1 ) 


Here N A - N A (t,T), is the number density of cosmic ray isotopes A = LiBeB at time 
t with energy (per nucleon) between T and T + dT. Note that eq. (1) is the general 
propagation expression for A a primary or secondary species, with different terms being 
important for each. The first term on the righthand side, J A , includes all sources for A: 


MT) = QAffl + Znj Q c dT'MT') d -$(T,r) ■ ( 2 ) 

Here Q A (T) is a possible galactic source of A, given as a number rate per volume and 
per unit energy; and spallation production of A appears in the sum runs over projectiles 
* and targets j. Additionally, <j> t = N iVi is the cosmic ray flux spectrum of species i and 

18 the CTOSS sectlon for the process i + j — * A H . The second term of propagation 

equation (1) energy losses to the ISM, with b A = -(dT/dt) A allows for ionization from 

the Galaxy. The third term of the propagation equation accounts for catastrophic cosmic 
ray losses, 

Te// -E^vi «an, + — (3) 

with crhf' encoding spallation losses of A in the ISM, and r esc (T) being the lifetime for 
cosmic rays against escape. 

The propagation equation (1) is solved for the case of a steady state, dN/dt = 0, 
m which cosmic ray production is in equilibrium with the losses. One thus obtains the 
spectrum N of these elements, propagated from their source J. At present, we will ignore 
losses due to inelastic nuclear collisions, valid when r'J > n t c rj^v A (but see section 6.1). 
We assume that the primary cosmic ray species p, a, and CNO, have some homogeneous 
galactic source J — Q, and negligible spallation production or losses: a - 0. Writing 
the solution for these in terms of the cosmic ray flux </> = Nv, we have 

1 f°° 

^ T ^ = MTjjT dT '^ T ') ex P(-l R i(T')-Ri(T)]/A) (4) 
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TrfL 9 *' = Qi/piSM ’ and *' = P’ CNO. Also, i^(2U) = 
Jo ai l\pi /OA ) lSM is the ionization range which characterizes the average amount 
of material a particle with energy T A can travel before ionization losses will stop it 
(expressed in gem 2 , as X = pisMvt). The ionization ranges were taken from North- 
cliffe and Schilling (1970) for low LiBeB energies (< 12 MeV/nucleon) at which partial 
charge screening effects are important, and from Janni (1982) for higher energies, using 
the scaling law 

R a (Z-T) = A/Z 2 R p (T). (5) 

Here T the kinetic energy per nucleon and R p is the proton range. 

In figure 1 we plot the proton and a fluxes <j>i calculated from eq. (4). To show the 
effect of energy losses on the propagated flux, we also plot A q { , the solution of eq. (1) 
for negligible energy losses (6 t - & 0). As we will discuss below in sections 5 and 6.1, 
these losses are important at low energies and negligible at high ones. This is manifest 
in figure 1, which shows the scaling (f> = Aq to be followed closely at high energies, while 
at low energies the ionization energy losses significantly reduce the propagated flux from 
this scaling. This behavior is qualitatively similar for the two source types we consider; 
propagation differences at low energies are discussed in section 6.1. 

In solving the propagation equation for the secondary elements LiBeB, the source 
term J is assumed to have no primary source component, i.e. <5 = 0. The exclusive 
production of these elements occurs via spallation between the primaries and the ISM: 
<Ti J ^ ^ ne can wr ^ e an expression like eq. (4) for LiBeB, giving their cosmic ray 

spectrum. What we wish to know, however, is not this steady-state spectrum but 
instead the amount of LiBeB thermalized and added to the ISM. To compute this from 
the LiBeB spectrum one assumes that all such nuclei below some threshold kinetic energy 
Ttherm are thermalized and then one examines the LiBeB current b A N A (T therm ) below 
this threshold. Below the lowest spallation threshold there is no source term in eq. (1) 
for LiBeB, and the ionization loss term is much larger than the escape term. Thus to a 
good approximation the propagation equation (1) reads 

Ofb A N A (T) = 0,T<TT (6) 

and so the LiBeB current b A N A (T tkerm ) is constant for T therm below spallation production 
thresholds. One may thus choose any T therm < T t n h uc at which to evaluate this current; 
we choose ours to be right at threshold. 

With this method of computing the production rate of LiBeB via evaluation of the 

subthreshold LiBeB current, we can write the rate of LiBeB accumulation in the ISM 
as 

— d r = T Vi(t) jT. JT MT, t) 4(T) S a [T a (TI t] (7) 

ij 1 th 

where A — Li, Be, B, y A = n A /nu , the <f>i come from eq. (4), and we have ignored 
the small time variation of n^. 
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. , 6 , fa , C . t , 0r Sa ( Ta ’ *) acc °unts for the energy loss of nucleus A in the ISM and gives 

the probability of its capture and thermalization; it is a function of the lab energy T A 
of the daughter nucleus A, as well as the epoch t, and is given by 

S a (T a , t) = exp [- {R a {T a ) - R A (T th „ m )} /A (i)] (8) 

These energy losses in the ISM serve to trap the spallation products and so to allow 
h “ ,° contrlbu te to the LiBeB finally thermalized in the ISM. This process competes 
wi h the the spallation products’ escape from the galaxy, quantified by the escape rate 

T ° ,c 0t (1) ; Th,s 9 uantity a PP ea ™ in eq. (8) through A = PlSM VT esc the (possibly) 
energy-dependent average pathlength 1 , also in g cm' 2 . Note that as r„ c sets the scale for 
the cosmic ray residence time before escape, A is the amount of matter traversed before 
escape, s any therm < T th gives R(T t h „ m ) <C A, we may put S A ~ exp(—R A / A). For 
forward kinematics, with a light cosmic ray nucleus impinging on a stationary, heavy 

„ ^nucleus, ?? IS sma11 for ty P lcal spallation energies, and so R A < A and S A ~ 1. 

or inverse kinematics, with a heavy cosmic ray nucleus on a light ISM particle, S A 
can differ significantly from unity and the LiBeB yields are reduced accordingly 
Following Walker et al. (1992), we put 


<H(T,t) = y? R {t)q p (T,t) 


( 9 ) 


with i - p a, CNO: i.e., we posit the constancy of the cosmic ray isotopic and elemental 
ratios at the source and over the entire energy spectrum. We then choose to make the 
more serious assumption that we may express the proton source strength in the separable 

9p(T,t) = q„(T) f(t). (10 ) 

We wdl in fact consider various Population II source spectra q (see section 6.2), but we 
will not allow a fully general energy or time dependence. While this has the immediate 
f, V “ age ° slrr| plifying the calculations, one may view this assumption as a postulation 
that the present mechanism of cosmic ray acceleration does- not differ dramatically over 
ime m its energetics, but only in its net cosmic ray output. 

In our analysis of cosmic ray model features, and in the accompanying figures, we 

will concern ourselves with the effect on the LiBeB production rates-or rather their 

ratios as calculated by eq. (7). Given a set of cosmic ray and ISM abundances, and 

a confinement parameter A, these rates may be evaluated numerically. As in previous 

empirical values of partial cross sections tabulated in by Read and Viola 

(1984). To proceed further and integrate the rates to get the LiBeB yields would require 

a lull model for cosmic ray and chemical evolution. We discuss this issue in the next 
section. 

conlLTn tit* rTa ^ eXP ° nential in 88 » eq. (8), we are assuming A to be 

constant in energy. For further discussion see section 6.1. 


6 



As a cosmic ray model of propagation, the above amounts to the simple leaky box 
(reviewed by, e.g, Ceasarsky 1980, 1987). This is the simplest model used to describe 
cosmic rays today, and as such has obvious computational advantages. While the sim- 
plicity of the leaky box makes it a useful tool, it is physically unrealistic and thus in 
some cases quite inaccurate. A proper treatment of cosmic ray propagation must ex- 
plicitly include diffusion effects only sketched with the leaky box confinement parameter 
A. Also, a proper model must abandon the simple leaky box assumption of spatial ho- 
mogeneity of sources and interstellar material. Despite these shortcomings, however we 
will follow previous authors and adopt this model in what follows, both for the above 
pragmatic reasons and moreover because the epoch we consider is so poorly understood 

in its relevant details that adoption of a more detailed description is not warranted at 
the present time. 


4 The Role of Chemical Evolution 

Although we can calculate the rates in eq. (7), for a given set of parameters, using 
well-understood physics, to integrate these rates requires knowledge of the chemical and 
dynamical evolution of the cosmic ray and ISM species in the Galaxy’s early history 
buch knowledge is sketchy at present, but some reasonable assumptions will prove fruit- 

If we assume that the H and He abundances do not change much from their values 
as set by the big bang, we have 


m 

VHe{t) 


H BB 

,.BB 

VHe 


0.08 . 


( 11 ) 


If we further assume that the CNO nuclides evolve at the same rate, so that their ratios 
remain constant within the early epoch considered here, we put 


Vc(t) _ VN(t) _ yo(t) 


y° BS 


,,055 

Vn 


9I obs 

Vo 


( 12 ) 


assumed to hold for times t less than the galactic age r at the birth of the Be or B star 
having CNO abundances denoted here as OBS. One should interpret equation 12 with 
some care, as approximation of constant C:N:0 ratios is not strictly true. We expect 
these elements to have different sources: O is made in type II supernovae but C and N 
are made primarily in intermediate mass stars. Thus any differences in the evolution of 
these sites will result in differences in the abundance ratios. Fortunately, large differences 

upma C r N and 0 production on *y important at early times, and the observed Pop 
11 CNO abundances are consistent with constant C:N:0 over the timescales important 
here. One should note, however, that the O/CNFe ratio, though the roughly constant, 
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does exceed the solar ratio. This comes about because of the predominance at this 
epoch of type II supernovae over type I’s. The effect is important for our purposes, as 
spallation off of O gives a lower B/Be ratio than spallation off of C. 

We again simplify by extrapolating the present relation between abundances of cos- 
mic ray and ISM species, namely we assume that the relative enhancement of the cosmic 
ray CNO abundances over interstellar CNO abundances has remain constant in time 


y? R (t) _ (y? R \ 


Vi(t) 


Vi ) 


present 


(13) 


where i - CNO and the righthand side denotes the present value. We will in fact adopt 
the stronger assumption, made in most work, that 


y cn (t) = y ISM (t) ( 14 ) 

but we note that this assertion is false for the present cosmic rays. In particular, Asaki- 
mor, et al (1993), and references therein, report that found a H/0 depletion of about a 
factor of 2 relative to solar. The proton-to-helium ratio in the cosmic rays is energy de- 
pendent and varies from ~ 0.2 at the lowest measured energies to ~ 0.05 at the highest 
energies which interest us here. Additionally, as shown in Buckeley et al (1993) He/O 
in the cosmic rays is depleted by a factor of ~ 4.6 relative to its local galactic value 
In this paper we have normalized to the protons, and thus adopting the observational 
results would amount to reducing y He by 2.3, and enhancing y CN0 by a factor of 2. We 
wi instead use eq. (14) to allow comparison of ours with previous results. Bear in mind 
however, that the cosmic ray abundance scalings are uncertain by at least a factor of 2 
One Proposed model which clearly contradicts this assumption is that which explains 
heBe/O constancy by having the CNO either as targets or projectiles be at supernova 
production abundances. We may model this behavior by putting 

y C o R {t) = y s o upernom ~ J £«— » yOBs (15) 

In what follows, however, wherever not explicitly noted otherwise, eq. (13) will be as- 
sumed to hold. 

If we allow the confinement pathlength to vary with time, then as we will see below 
we may not further simplify the rate equation (7), or the integral that is its solution.’ 
the confinement is taken to be constant in time— a dubious proposition, as we shall 
see— we may now separate the propagated flux: 


^ . ( 1 6) 
This afiows one to write the solution to the rate equation as a sum of terms, each of 

which has a factor that involves an integral over time and a factor involving an integral 
over energy: ° 6 

y A (r) = £ Aij ($<xd) (17) 
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with the “exposure time” given by 

A i; (r) = ^ dt V j(t) yf R (t) f(t) (18) 

with r is the age of the Galaxy at the birth of the star in question, and the “reduced 
rate is 

A f°° 

= J Tnuc dT UTy t) 4{T) Sa(T a , t) . ( 19 ) 

th 

Note that these factors are different for forward and inverse kinematics, that is, ^ 

\® a ij)r- As mentioned above, this stems from the longer range and greater chance of 
escape for the fast A-nuclei produced in inverse kinematic collisions. 

With the assumptions made thus far in this section — most importantly, that of the 

^ of A ’ the number of independent exposure times may be reduced to three 

(bw, WSSOF): 

f A a>a % ~~ j Oi 

Aij = < A ^q 1 forward kinematics (20) 

1 a$ inverse kinematics 

Of these, we will take A 0 , a as the fiducial quantity, as it measures the integrated flux 
enhancement exclusively (ie with no explicit dependence on chemical evolution) 

= Jo dt f(t) (21) 

We then define the quantity {O/H), which is a measure of the average abundance of 
oxygen relative to hydrogen,: 


{O/H) ism 


A 

fo dty^)f(t) 

lo dt f(t) 


( 22 ) 


with a similar expression for (0/H) OR . Note that (O/H) is to the weighted average 

ot O/H at time r. This average weighted by the net flux enhancement fit), and is 
normalized. 

With this notation we have 


yB, 


•,b A 00 I (0/H) IS m 52' ($<r?j e,B )/ + ( 0/H)cr 52'(^ tr § e ’ B )r 

\ i) ij 


Vu = A< 


<*0+ i (o/H) ISM 52' (*<%)/ + (o/h) or £'<*<#>, 


(23) 


(24) 
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where the prime on the summation indicates that any i=j = a term has been deleted. 
Note that the ratio of LiBeB production rates, i.e. of eqs. (23) and (24), does not depend 
explicitly on A aa . However, the flux enhancement appears via the strong dependence 
on (O/H), which includes a factor of A aa in its definition. 

From eq. (17) it is clear that chemical evolution effects are encoded into the LiBeB 
abundances in terms of exposure time, in which these effects are convolved with cosmic 
ray flux evolution. While this means that the absolute abundances are dependent on 
chemical evolution, it also means that the relative abundances are much less so if the 
LiBeB isotopes compared originate in similar processes. This is the case with Be and B, 
which are both produced in pure spallation events, as opposed to Li, which is produced 
not only by spallation but also in a fusion events. This point was stressed by WSSOF, 
who noted that the ratio B/Be is the least dependent on our knowledge of the early galaxy 
and of early cosmic rays. WSSOF showed that any model for early cosmic ray LiBeB 
production yields a B/Be ratio lying within the range 7.5 < B/Be < 17. One might 
compare this with the average of the observed values in the most extreme Population 
II dwarf, HD 140283, which has B/Be = 6 db 3. This ratio is the strongest prediction 
of the cosmic ray model. However, while the theoretical prediction is firm, particularly 
in setting a lower bound to the B/Be ratio from cosmic rays, the experimental ratio 
is still difficult to obtain reliably. As discussed in section 2, the true errors to the 
measured value, when including systematic uncertainties, are certainly larger than the 
simple weighted mean we quote here. Comparisons with the observations are thus not 
all that compelling at the present time. 

Under the assumption that there is no additional (primary) source for 11 B, the B/Be 
ratio is weakly dependent on chemical and cosmic ray evolution. However, it is well 
known that the predicted value of the 10 B to n B ratio by GCR nucleosynthesis yields a 
ratio of about two, whereas the observed ratio in meteorites is just over four. Dearborn 
et al. (1988) proposed a possible thermonuclear production of n B in Type II supernovae, 
and Woosley et al. (1992) had suggested that neutrino processes during supernovae may 
yield a significant abundance of n B. Both suggestions are tested in a GCR and chemical 
evolution model by Olive et al. (1993), where it was found that while the “B-^B ratio, 
could be brought into agreement with the observation, the supernova sources, added in 
a significant primary source for U B which upsets the constancy of the B to Be ratio as 
a function of metallicity. Though the data seem to disfavor this primary source which 
predicts a B/Be ratio in excess of 50 at [Fe/H] < —3, the scant amount of data at present 
can not conclusively eliminate this possibility. 

In contrast, the ratios of Li with Be and B are always strongly dependent upon 
chemical and cosmic ray evolution because of the a — a production channel for Li. In 
addition, the dominant component of the observed 7 Li is primordial, further complicating 
comparison of Pop II LiBeB with GCR nucleosynthesis predictions. 
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5 The Be and B versus Fe Relations 

Many have suggested that since Be and B arise from the interaction of stellar nucleosyn- 
thesis products with cosmic rays, they should vary quadratically with elements, e.g. 0 
which form their targets. This dependence is argued to follow from the fact that cos- 
mic ray sources are presumed to be supernovae which are the producers of the oxygen 
targets. Of course such an assumption is extremely naive since cosmic rays have a large 
confinement time and hence the instantaneous recycling approximation that yields this 
quadratic dependence is quite inappropriate. Even more significant is the fact that the 
data does not show this dependence and is instead best fit by a Be slope relative to O 
t at is much closer to 1 than to 2. This discrepancy has received some attention and 
has even been mis- interpreted it as a failure of the cosmic ray model. It is obvious that 
reports of the death of early cosmic ray LiBeB are, however, greatly exaggerated 

Let us now go through the logic in greater detail. The argument for a quadratic 
scaling of Be and B with 0 proceeds from these assumptions: (a) Type II supernovae are 
the dominant source for 0, so one expects a relation to N sn , the total integrated number 
of supernovae, of the form 0 oc N sn . ( b ) Supernovae are also likely to be the acceleration 
site for cosmic rays, so one might naively expect nthe cosmic ray flux (f> to vary with 
the current rate at which these acceleration sites are created: <f> oc dN SN /dt if we 
ignore confinement and assume instantaneous recycling of the cosmic rays. Put slightly 
differently, in the leaky box model we have adopted, we have assumed an equilibrium 
to hold between the rate at which cosmic rays are lost and the rate at which they are 
created. Thus in this model one can argue that the cosmic ray flux is indeed proportional 

to this creation rate. Then, taking the essential pieces of the rate equation (7), we obtain 
the naive result 

dBe/dt oc 0 <f> 

oc Nswd Nsn / dt 

=* Be « ° 2 (25) 

which gives the advertised quadratic dependence, with an analogous equation for B 
evolution. 

While assumption (a) seems solid, assumption ( b ) is not because it accounts for 
evolution of cosmic ray sources but not for changes in confinement. Indeed the evolution 
of A is the critical feature allowing one to avoid the naive secondary-to-primary ratio 
for the B- or Be-to-0 ratio of eq. 25. For a given source density q of primary cosmic 
ray elements (at some instant f), the propagated flux <j> increases with A, as at larger 
A the cosmic rays traverse more matter before escape and so the steady state between 
production and escape is reached at a higher flux. In the limit of no energy losses, 
valid at high energy, the relation is <f, = A q. With energy losses, this relation is energy 
dependent, but for the total flux PCV find that it varies roughly as ^ oc A 1 / 2 . Thus we 
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should alter the cosmic ray flux scaling to be <f>(t) ex A' /2 (t)N SN (t), and it follows that 
changes in the early confinement alters the cosmic ray flux, and its associated LiBeB 
production. Specifically, in the place of eq. 25, we have the more general relation 

dBe/dt <x A n OdO/dt 

=>dBe/dO oc A n O (26) 

with nle 1 the index for the scaling law between cosmic ray flux and confinement. If A 
can vary with time — and indeed it is hard to imagine how to avoid such variance— then 
the simple quadratic relation between Be or B and 0 is clearly not the prediction of 

the cosmic ray model. The effect of the variance of A is further discussed in the next 
section. 


6 Cosmic Ray Model Dependences 

We wish to investigate the dependence of calculated elemental ratios Li/Be and B/Be 
on various model parameters, for different epochs, i.e. different [Fe/H]. Such a calcula- 
tion, however, requires specific models for galactic chemical and cosmic ray evolution. 
Without constructing such models, what we can examine, at a given [Fe/H], are the 
ratios of production rates dLi/dBe and dB/dBe and their dependence on escape length. 

This comparison of production rates provides useful bounds on the elemental ratios 
any full model would predict. Specifically, at a given epoch [Fe/H] and confinement 
A, dLi/dBe gives the value Li/Be would take if the ISM abundances and cosmic ray 
confinement had been constant for the entire history of the galaxy up to that epoch. 
The ISM abundances of heavy elements CNO certainly are not constant but build up 
from zero, while the ISM abundances of H and He do remain roughly constant. The 
chemical evolution of the ISM affects Li and Be or B differently; Li is largely unaffected, 
being made mostly by a + a fusion, while Be and B are made by spallation of the 
CNO targets whose abundances will have always been rising to the levels observed at a 
particular [Fe/H]. The calculated dLi/dBe will thus provide a lower bound on the actual 
Li/Be a full chemical evolution model would predict for a given [Fe/H]. 2 The B/Be 
ratio, however, does not change at all with CNO evolution as long as the C:N:0 ratios 
remain constant, and thus is equal to dB/dBe . These ratios are roughly constant over 
the Pop II epoch we consider here, but they do of course change eventually as the 0 
overabundance diminishes to it Pop I value. 

2 Just as the CNO abundances will have changed, so too the cosmic ray flux strength and confinement 
will have evolved. Though an overall change in flux strength does not affect elemental ratios, evolution 
m confinement will change the spectrum, with high early confinement leading to a harder spectrum in 
the past (see next section). This decreases the Li/Be ratio, as the a -f a production falls off at high 
energies. This goes in the opposite direction as the chemical evolution effect. However, we will in 
general explicitly examine the allowed confinement parameters, and in this context we may'understand 
the dLi/dBe ratio to be bounded from below. 
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6.1 Confinement 

p CV h av e detailed a chemical evolution model that can fit observed Be and B data 
and does produce a slope versus iron that is smaller than the naive quadratic depen- 

ence.^ In their model, the pathlength is assumed to be larger than the present ~ 10 
g cm due to the larger confinement volume of the halo phase versus the disk phase 
as well as the larger early gas fraction. They suggested that the confinement varies 
over the range 10 gem" 2 < A < 1000 gc m - 2 . This variation arises specifically from 
scaling the escape length as A <x <r „H; where the gas mass fraction <r g decreases with 
time according to the adopted chemical evolution, while the disk scale height H has an 

exponentially decreasing evolution imposed on it by hand, to simulate the dynamics of 
the halo collapse. 

While this line of reasoning is attractive, one must be cautiously aware of the uncer- 
tainties involved. For example, the early behavior of the main confinement mechanism, 
the galactic magnetic field, is not well understood. The favored means for field gen- 
eration is through dynamo action tied to the galactic rotation, setting a timescale for 
generation of r rot ~ 10 8 yr. Note that in this case the galaxy does not as efficiently 
retain LiBeB long enough to thermalize them, particularly given the flatter (harder) 
cosmic ray spectra and/or the inverse kinematics some have suggested (see section 6.2. 

Furthermore, Malaney and Butler (1993) have pointed out that the escape length is 
bounded from above. Physically, this happens when the confinement becomes so large 
that nuclear destruction of fast LiBeB becomes important. In terms of eq. (1), we have 

I- 1 1 

A - Aon/ + AZ W 

where A = £i n^ 1 / p !SM . Malaney and Butler note that A„„„ is species-dependent 
(whereas A cmf ,s not), and they estimate A nuc < 200 gem' 2 for protons, thus setting 
an upper limit for A. Prantzos and Casse (1993) have included this effect in the PCV 
model and find it makes little change in their results. 

What is the affect of a changing confinement? A variation in confinement, acts 
not only to alter the flux intensity but also to change the shape of the spectrum This 
behavior is apparent in figures 2(a) and ( b ), which plot spectra for different values of con- 
finement. Notice that in both plots the total flux is higher for larger confinement, but by 
different amounts in different energy regimes. The lowest energy (< 100 MeV/nucleon) 
portion does not enjoy the same enhancement as the high energy (> 1 GeV/nucleon) 
region. One can see in eq. (4) that the spectrum turns over roughly at the energy at 
which R(T) - A. A spectral peak at a larger T means increased power at higer ener- 
gies. In other words, increasing confinement enhances the highest energy regime and so 
hardens the spectrum. 

Note the different low-energy flux behaviors for the two source types, the forms of 
which are given explicitly below (section 6.2, eqs. (31) and (32)). The momentum source 
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diverges as T -* 0, whereas the energy source is finite for small T. These differences 
manifest themselves in the propagated fluxes. The flux derived from the momentum 
source has a mild turnover at small energies, and the different A fluxes all converge. 
Th!s behavior is contained in the flux equation (4). At low energy (T < m p ), the escape 
ength R(T) ~ T/ (dT/dX) < A, is slowly varying, and so the exponential in eq. (4) 
changes slowly as well. On the other hand, the source is a featureless and steeply falling 

power aw. q(p) ~ p ~ T . Thus one can, at sufficiently low energy, ignore the 
confinement parameter altogether and obtain 


1 r°° 

dT '^ ( 28 ) 

which is independent of confinement. A source power law in total energy, however, 
the source has a characteristic scale m p and does not strongly vary over this scale.’ 
bince K(1 )/A does change appreciably in this range, eq. 28 is not appropriate, and 
one must use the full flux equation, with its the strong dependence on confinement. 

hus the propagated fluxes do not display the same convergence as those derived from 
a momentum-law source. 

The well-known (e.g. PCV) dependence of the propagated flux on confinement has 
important consequences for elemental ratios. Although changes in the total cosmic ray 
intensity will not affect elemental production ratios, spectral changes can. One calculates 
LiBetf production rates by convolving cross section with flux. Spallation and fusion cross 
sections genencally have low energy resonance peaks up to a few hundred MeV/nucleon 
there is a high energy asymptotic behavior characteristic of the interaction. Namely’ 
the cross sections approach a constant at high energies for spallation interactions with a 
many body final state, but suffer an exponential dropoff for fusion reactions which have 
ew-body final states. (It should be cautioned here that the continued exponential falloff 
of the a + a fusion reactions above a few hundred MeV is only an extrapolation since no 
data exists at higher energies; c.f. section 6.3.) Thus the effect on dB/dBe —the ratio 
of pure spallation products— is to change the emphasis on the resonance region relative 
to the plateau. Both p and a cross sections for B production have lower thresholds than 
those for Be, and B cross sections also display higher peaks relative to their plateaus. 
Consequently, the dB/dBe ratio is lowest at cosmic ray energies in the plateau region. To 
produce a low B/Be ratio requires a spectrum that gives a larger weight to this energy 
regime, 1 e. a harder spectrum. Just such a spectrum arises from a large confinement 
scale. The absolute effect on the dB/dBe is small, but important given that B/Be is 
constrained theoretically and observationally to lie in a small range. 

For the dLi/dBe , B production ratio, spectral flattening via large confinement leads 
to more pronounced changes. The effect again stems from the high energy cross section 
behavior: while the Be and B cross sections are asymptotically flat, the a+a Li fusion 
cross sections are assumed to fall off exponentially at high energy. Thus harder spectra 
have a strong effect dLi/dBe as they emphasize an energy regime where the Li production 
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is not just reduced relative to Be, but in fact negligible. Thus dLi/<ffie strongly decreases 
with increasing confinement. 

We have computed numerically the ratios of rates for the range 10 gem -2 < A < 
1000 g cm 2 . The results appear in figures 3 and 4, and are in agreement with our expec- 
tations. The strong dependence of these ratios on the different forms of source spectra 
will be further discussed below (section 6.2). In reading these figures and comparing 
to observational data, one must bear in mind the need to include evolutionary effects 
properly in order to go from the plotted ratios of rates to elemental ratios. Furthermore, 
as argued above, we see the falloff of both dLi/dBe and dB/dBe at large A. Note that 
the scales in these figures are very different. As A is increased from 10 to 1000 gem -2 , 
dLi/dBe falls off by a factor of 10, whereas dB/dBe drops off by only about 15% for the 
momentum source spectrum and even less so for the energy source spectrum. 

In figure 5 we plot the production ratio dLi/dBe against [Fe/H]. Such a plot sug- 
gests a comparison with observational data, but bearing in mind the caveat that the 

dLi/dBe value is only a lower limit to the actual value a full chemical evolution model 
would give. 

In our considerations thus far we have followed most previous work on early LiBeB 
in putting A to be a constant in energy. This parameterization offers a simple way of 
trying to characterize the early galactic conditions which, as has been pointed out, are 
poorly understood. While offering simplicity, this approach does suffer the criticism 
that the choice of a constant escape length does not fit the observed present day cosmic 
ray pathlength distribution. As shown most recently in Garcia-Munoz et al. (1987), the 
escape length varies with energy, increasing slowly to a maximum around a GeV/nucleon, 
then decreasing as T °' 6 . They find that their empirical form for A(T) is well fit by 


A(T) = A 0 


(T + rn^y 
\7i -f m P J 


1 — 0.2 exp 




(29) 


with A 0 = 11 gem 2 , T, = 850 MeV, and T 2 = 300 MeV; /? = +2 for T < T u and 
/? = -0.6 for T > 1\. A plot of this function appears in fig. 6. As is well-known, an 
energy-dependent pathlength changes the solution to the propagation equations (1 and 
7), with the replacements 


R ( T )/& + f Q w ( T ,j^ T ,j (30) 

where w(T) = (dT/dt)/p ISM v. 

We have allowed the normalization parameter Ao in A (T) to vary, and figures 7 
and 8 compare the resulting dLi/dBe and dB/dBe to the same ratios obtained with 
a constant A (T) = A 0 . Note that the changes are relatively small. This is because 
the energy dependence of the pathlength is most important at high energies, above the 
region R(T) ~ A(T) at which the flux turns over. However, the majority of the flux 
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resides in the region around and below this maximum. Since these energy regimes are 
largely unaffected, the LiBeB production is weakly affected as well. For the momentum 
source spectrum the results with and without an energy- dependent pathlength differ by 
a constant factor of roughly 25%, while for the energy source spectrum the difference 
varies but is always small < 5%. The smallness of the difference between the two 
confinement behaviors justifies the A(T) = A 0 approximation used in this work and 
elsewhere. The systematic decrease of the yields for a given scale A 0 follows from the 
energy dependence, under which A has a maximum, at T = T u of only 0.8A o . Thus the 

confinement is always less than this value, and very much lower for the important low 
energies. 

We have investigated the effect of inelastic nuclear reactions between the primary 
cosmic rays and the ISM on the effective confinement given in eq. (27). As discussed by 
Malaney and Butler (1993), this effect contributes an extra species-dependent and energy 
dependent term to the escape length. Because the pathlength is energy dependent, one 
must use eq. (30) to compute its effect. We have done so, using cross sections given in 
^i 72 ^’ Townsend & Wilson (1985), and Bystricky, Lechanoine-Leluc, & Lehar 
(1987). The results appear in Figure 9, and the behavior is as expected. In the usual 
case, dLi/dBe decreases with increasing A e// = A eac , and we now have A e// < A e5C . 
Thus, we expect dU/dBe to be consistently larger at a given A esc , as is the case. This 
effect is more pronounced at high A esc , where the nuclear loss term dominates. 

Another property of the early galaxy with bearing on confinement is the state of 
ionization of the ISM. Energy losses are increased in an ionized medium, leading to a 
shorter ionization range. Because the confinement always appears as Rf A (for A(T) 
constant in energy), this effectively increases the escape length PCV account for ion- 
ization by doubling the energy losses for neutral media, a procedure which amounts to 
assuming a fully ionized medium (Ginzburg & Syrovatskii 1964). If indeed the ionization 
state evolves during the Pop II epoch, the effective change in A can affect the cosmic 
ray spectrum independent of the considerations in eqns. (27) and (30). 

6.2 Spectrum 

Both WSSOF and PCV have independently argued for a flatter spectrum in early cosmic 
rays than that presently observed. The key observation supporting this argument is the 
B/Be ratio. Although the observational uncertainty here is clearly large, the current 
preliminary numbers are low compared to the present day observed value of 15 ±3 (with 
e.g. Mathews, Walker, and Viola (1985) predicting ~ 17 for Pop I), and to the Pop II 
value of ~ 14 calculated in WS (using a transported spectrum 6(T) oc (T + m )~ 2J ) 
and Pop I C:N:0 ratios. ^ p> h 

The argument for a flat early cosmic ray spectrum arises if one takes a low central 
value for the Pop II B/Be ratio seriously. (Remember values below 7 are impossible 
m an y spallation model, and so far hypothetical stellar processes enhance B.) Be and 
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B arise from the same spallation processes, so there are a limited number of available 
differences to exploit in fixing the ratio. Their differences in mass and charge lead to 
different ionization losses and ranges, and consequently different confinement for the two 
elements (at a given energy per nucleon). However, with the ionization ranges scaling 
a/w of eq (5), it is clear that R B (T ) < R Be (T ), and thus energy loss and thermalization 
effects (c.f. eq. 8) favor B rather than Be and so act to increase the B/Be ratio. The 
other important difference between Be and B production rates lies in their spallation 
cross sections. As argued in section 6.1, the differences make the B/Be ratio sensitive 
to the cosmic ray spectrum, with a lower B/Be for flatter spectra. 

Another argument for a flatter spectrum has been offered by PCV (but weakened by 
Malaney and Butler’s (1993) arguments; see Prantzos k Casse (1993)). PCV note that 
the low observed Li/Be ratio (once the primordial Li component has been subtracted) is 

more easily reproduced by a flat spectrum because of exponential dropoff of a -f- a cross 
sections. 

In addition to the effect of changes in the confinement, the propagated spectrum 
obviously depends on the source spectrum. PCV choose as a source spectrum a power 
aw m momentum, q(p) oc p 2 , a type of spectrum appearing in models of cosmic ray 
shock acceleration. Such models assume acceleration via scattering off of supernova 
shock waves, and have been successful in reproducing galactic and heileospheric spectra 
m detail (c.f. the Blanford k Eichler (1987) review). Note that the PCV spectral index 
of 2 is the universal limiting value for very large Mach number shocks. We note, however, 
that despite the success of and active interest in this model, there remains disagreement 
upon the the appropriate model for acceleration. It is thus worthwhile to examine the 
sensitivity of the LiBeB yields to the assumed spectrum, where there remains some 
flexibility to choose the spectral index. 

Most important for LiBeB production is the choice of spectral type. While shock 
acceleration models suggest source power laws in magnetic rigidity and hence momentum 
are good candidates, the measured (and hence propagated) cosmic ray flux at earth 
is also consistent with a source spectrum in total energy per nucleon. The two are 
of course difficult to distinguish at the relativistic energies at which the cosmic ray 
measurements are unaffected by solar modulation effects. Furthermore, solar effects 
introduce large uncertainties at lower energies, rendering unreliable the data that could 
make this distinction. In addition, a power law in total energy has the advantage that it 
remains finite at low energies, thus providing a mockup of physical processes responsible 
for preventing a divergence in the source at low energies. 

Thus we have tried source spectra of power laws both in momentum and in total 
energy per nucleon, i.e. 

q(p) (X p~~' 

=* <?( T) <x(T + m p )(2m p T + T 2 )-< 1+1 >/ 2 , and 
q(T) oc (T + m p )^ 
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(31) 

(32) 




with m p the proton mass. For the shock acceleration model-inspired momentum spec- 
trum, the index 7 is given by 

7 “( 2 + ife) W 

for the case of large Mach number M = v sho ckl c sound , the ratio of shock velocity to the 
sound speed in the shocked gas. Note that the Mach number depends upon the local 
ISM sound speed, which in turn depends upon, e.g. the ISM density and pressure. If 
the early galaxy, or its primordial protogalactic building blocks, were to be significantly 

hotter or denser than today, then potentially the cosmic ray spectrum could be quite 
different. 

For the total energy spectrum we take 7 = 2.7 , which is consistent with the mea- 
sured, transported cosmic ray spectrum. Note however, that this choice is inappropriate 
for the simple leaky box with energy dependent escape length. In this case, we live in 
t e confinement region (the galaxy), and so the escape energy dependence is folded into 
that of the source. The proper source index is thus y, ource = %ts - ,3 cs 2.1, with 6 the 
high-energy index of the energy-dependent A from eq. (29). This is not the case in, for 
example, the nested leaky box model, in which we do not live in the confinement region 
(which shrouds the sources), and so we do not see its energy dependence: 7 source = Jobs • 
Again, the model uncertainty argues for examination of a range of source indices. 

Note the dramatic effect on the Li/Be and B/Be plots (figs. 3 and 4). The large low 
energy flux of the momentum spectrum makes this spectrum softer than one in total 
energy, and the resulting elemental ratios behave as one would expect for such a soft 
spectrum. As discussed above in section 4, low energy resonance behavior is emphasized 
and leads to large Li/Be and B/Be ratios. 

We have allowed the spectral index 7 to vary over the range 2 < 7 < 3. Results 
appear in figures 10 and 11 from which it is apparent that the LiBeB production rates 
from a source with a momentum spectrum are more sensitive to changes in spectral 
index than are the rates for a source total energy law. Indeed it is clear from fig. 10 
that for source spectra in momentum, indices with 7 > 2.5 lead to Li-to-Be relative 
production rates in ratios approaching 10 3 . This is clearly undesirable given that all 
observed Be abundances are > 10 -13 , thus implying associated Li production larger 
than the observed plateau, even without the large primordial Li component. 

While the relative B and Be rates (figure 11) do not vary as dramatically for the 
momentum spectrum, here too they are much more sensitive than the same ratios given 
a total energy spectrum. Additionally, the same range of indices that overproduce Li 
relative to Be give the largest values for B relative to Be as compared to low values for 
the current preliminary but uncertain data. Nevertheless, it is intriguing that the range 
of indices that fail for the Li to Be ratio seem to also do the worst for B to Be. Thus 
momentum spectra with indices steeper than about 2.5 are less favored at the present 
time. Note that we get no such constraints for source spectra in total energy, as we see 
only mild sensitivity to a source index for Li to Be (for which all values are completely 
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acceptable), and completely insensitivity for B to Be. 

In figure 12 we plot dLi/dBe vs A, as in figure 3, but with different spectral indices, 
bor each spectral type, the production ratio grows with the steepness of the spectrum, the 
growth being more pronounced for the scale-free momentum law. Figure 13 is a similar 
p ot or the fiB/dBe ratio. Note that this ratio is mildly sensitive to the momentum 
spectral index, which can change dB/dBe by 1 to 2 units, i.e. 10-20%. Such a change 
rrnght be detectable as observations improve. The dB/dBe ratio is even less sensitive to 
the index of a source with a total energy spectrum. The change over the whole range is 

As we have alluded to earlier, a direct comparison of the data and the model cal- 
culations given here is extremely difficult but can be suggestive. The ratio of of Li! Be 
given in table 1 contains the primordial Li component as well as the cosmic-ray produced 
Li. In some cases this correction might be quite large (a factor of 10 or more) so that 
the ratios in the table represent extreme upper limits. If we had a specific cosmic-ray 
nucleosynthesis model, one could, as shown in Olive and Schramm (1992) extract a Big 
Bang Li abundance. For the cosmic-ray nucleosynthesis model used in WSSOF this led 
to a value [Li] bb = 2.01 ± 0.07, however as we can see from figures 5a and 56, there 
is significant model dependence in the Li/Be ratio. Also, the absolute [Li] values are 
plagued by the poorly given systematic errors mentioned in sectios 2. Furthermore, as 
we have also noted, the element ratio will also be affected by chemical evolution. Thus 
one can with certainty only require the data points as given in the table to lie above any 
model curve in the figures. The trend in the data is also seen in the figures, ie. at high 
metallicities the Li/Be ratio is smaller. Note that a large primordial subtraction would 
avor the lower ratios found with the cosmic-ray energy source spectrum. 

To get a feeling for the comparison between data and these models, we can make a 
primordial subtraction corresponding to [U] BB = 2.0. It hard to make the primordial 
Li abundance much smaller (Walker et al 1993; Krauss k Romanelli, 1990) and if it 
were much larger there would be no room for cosmic-ray produced Li. In Figure 14 we 
show the values of Li/ Be with the above primordial subtraction. Error bars are too 
large to display. We also plot the ratio from the energy and momentum-spectra for 
A — lOOg cm as a function of [Fe]. 

6.3 Other Model Features 

The conventions for reporting spallation/fusion cross sectional data are strictly speaking 
incomplete for calculations of LiBeB production, and thus require a model to make some 
assumption regarding the kinematics of these interactions. Specifically, the tabulated 
data and semi-empirical fits give the total cross section for LiBeB production at a given 
incident energy. As seen explicitly in eq. (2), however, we require the full, differential 
cross section dcr (j - / dT Li BtB {TuBeB,T CR ). Meneguzzi et al. (1971) and many subsequent 
authors address this issue by assuming that the differential cross section is very sharply 
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peaked around those values of the LiBeB energy for which TLiBeB — T C r. They put 


dot ! 


(ITt b d ( TubcBiTcr ) — <yfj(Tcn) 8(T LiB eB — T C r) . 

LiBeB 


(34) 


with crfj(TcR) the tabulated cross section. 

While the approximation of sharp peaking over a small range of product energies 
appears to be borne out experimentally in the cases where it has been checked, the 
other approximation, that of equal energies (per nucleon) in the initial and final state 
is only valid for the case of heavy cosmic rays on light targets. We have relaxed the 
latter assumption, allowing for the CR energy to be related to the LiBeB energy via 
some function E/(Tcr), changing eq. (34) to 


derf- 


dT t 


LiBeB 


(TLiBeB^TcR) — vfj(TcR) S {T LiBeB ~ Ej(Tcr)) 


(35) 


With this more general formalism we have investigated two possible kinematic be- 
haviors. One ansatz is that in the center of momentum frame the LiBeB is produced 
at rest, with the lighter debris moving in such a way as to conserve energy. This is the 
assumption we have made throughout this paper. One should compare this to another 
common assumption that in the frame in which the light target is at rest, the kinetic 
energy per nucleon of the daughter nucleus is the same as that of the heavy CR parent. 
We have plotted the dLi/dBe ratio for these two cases in figure 15. 

The issue of spallation kinematics is raised if one considers an alternative to the 
usual production scenario we have presented so far. Duncan et al. (1992) and PCV both 
attempt to address the question of the Be vs Fe slope with the suggestion that perhaps 
most of the LiBeB production is not from light (p, a) cosmic rays on heavy (CNO) ISM 
nuclei, but the reverse: heavy cosmic rays on light ISM nuclei. Furthermore, if the heavy 
cosmic rays do not have the ISM abundance or something near to it, then this could 
yield significantly different Be and B production. Increased cosmic ray heavies would 
result, for example, if supernovae were to accelerate their own ejecta. The evidence is not 
definitive on this issue, but there are possible problems in maintaining consistency with 
an assumed first ionization potential dependence of the cosmic ray source abundances. 

At any rate, we have investigated the dependence of the spallation yields on forward, 
reverse kinematics. Figure 16 compares the dhi/dBe ratio for the cases of normal 
kinematics (i.e. LiBeB synthesis primarily from light CR on heavy ISM with a small 
admixture of the reverse) to reverse kinematics, in which we assume y§NO ~ Vcnoi 
and thus dominate the LiBeB production. Note that here the dLi/dBe ratio is much 
smaller, as the a + a process is not as dominant as in the usual early galactic case. 

Note that, in considering the thermalization of LiBeB, authors from Meneguzzi et al 
(1971) to PCV apparently use S = 1 over all energies in the case of light cosmic rays on 
heavy targets. These authors note that the daughter nuclei move more slowly than the 
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cosmic ray parent and so are likely to be stopped. However, as we are considering cases 
where the flux can be very hard, then we may still have quite energetic daughter nuclei 
and it is important to properly allow for incomplete stopping. Note that Li, which is 
made primarily by a + <* fusion, is thus the only spallogenic element we consider for 
which (in the dominant production process) the cosmic ray projectile not much lighter 
than the ISM target. Consequently, the Li will be the fastest daughter nuclei, and thus 
accurate treatment of trapping is particularly important. To compare with previous 
work, we have computed yields with 5 s 1 for both kinematics. Results appear in figure 

We see in figure 17 that the direction of the effect depends upon the spectral type. 
For a momentum source, one gets a higher dLi/dBe ratio when allowing for realistic 
trapping of LiBeB. The change, ranging from about 40% at large A down to 10% at 
small ones, arises because Be production at high energy exceeds that of Li (because of 
the a - fa falloff at high energy), and these high energy Be are less likely to be thermalized 
than slower daughters. This is not the whole story for a source law in energy. Here, 
the spectral peak is much more pronounced and at a high energy (~ m p ) than the 
(lower energy) peak of the flux from the momentum source. The flux from the energy 
source thus rises steeply at the (~ 100 MeV) dropoff in the Li production cross sections. 
Consequently, the majority of Li production occurs in the beginning of the high energy 
tail of the cross section. By disfavoring the high energy products through incomplete 
stopping, one removes the faster Li made at the highest overlap of cross section and flux. 

With 5—1, however, these fast daughters are included and so raise the Li production 
relative to Be. 

Let us now look at the dependence of the dLi/dBe production rate on the asymptotic 
value of cross section for a+a — > 6,7 Li. This is a fusion process and is physically different 
from the other spallation processes. The data shows that these cross sections decrease 
exponentially with energy, although the highest energy at which a definite cross section 
has been reported is only about 50 Mev/nucl. Evidence for the continued dropoff comes 
instead from upper bounds to the cross sections, which go out to ~ 200 MeV/nucl, with 
no plateau detected. While Read and Viola (1984) recognize this, they suggest in their 
cross section tabulation that a plateau might nevertheless eventually be found to exist, 
and some work (notably PCV) has employed such a feature, although the data does not 
justify such an assumption. We have examined the Li/Be ratio with and without such 
a plateau, and find that it can have a significant effect, particularly for cases of large 
confinement or spectral index near 2, as these emphasize the high energy events in which 
the differences are important. Figure 18 compares the dLi/dBe ratio with and without 
the cross section plateau. While the differences are not large for a momentum source 
(< 10%), for a source power law in total energy, the effect is about a factor of 2. Since 
we see no evidence as yet for a plateau, in all our other computations none is assumed to 
exist. However, as the effect of such a plateau is significant, we feel it is very important 
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for experiments to investigate the behavior of the a + a -> 6 - 7 Li at high energy. 

Finally, we note the important but not decisive contribution to the overall uncer- 
tainty of experimental uncertainties in cross sections and in energy losses. Cross section 
uncertainties are typically at a level of Scr/a ~ 10 - 20%. While these errors are large 
enough to be significant, they clearly are smaller than the other model effects we have 
noted above, consequently making the model choices important ones. Energy losses are 
better known, with Sb/b ~ SR/R ~ 5% at low energies (where they are unimportant), 

and < 2% at higher energies. These errors do not significantly contribute to our overall 
uncertainty. 


7 Gamma Ray Production 

The same passage of cosmic rays through the ISM that produces LiBeB must also create 
7-rays, as noted by Silk and Schramm (1993), Fields, Schramm and Truran (1993), 
and Prantzos and Casse (1993). The primary 7 -ray source from the nuclear cosmic ray 
component arises from inelastic scattering of cosmic rays on the ISM, with p + p -+ 
7 r 2 7 . Thus the 7-ray spectral peak is around m r /2. 

The total (isotropic) 7-ray production rate by the nuclear component of the cosmic 
rays, per second per cubic centimeter, is given by 

J fOO 

' Tth dT MT)tf p (T) (36) 

in an obvious notation. Note that the calculation is similar to that for LiBeB produc- 
tion. Indeed this process occurs under very similar conditions to Li production, as both 
processes derive from reactions between protons and a nuclei, which in the early galaxy 
are undepleted both in the cosmic rays and in the ISM. Thus we may expect to scale 
the 7-ray production to that of Li via <r{a + a - Li) to the <r(p + p - 7) cross section 
ratio. One should note, however, that such a simple scaling ignores the different energy 
dependences of these processes, and ignores the issue of the incomplete trapping of the 
Li fusion products. 

Because 7-ray production so closely follows that of Li, the model dependences will 
be similar. In particular, the effects of cosmic ray evolution will be similar, as will 
the uncertainties in scaling the 7 -ray production to the Be. In a subsequent paper on 
specific galactic evolution models we intend to quantify the 7 -fluxes and use them as a 
complementary probe to the LiBeB abundance behaviors described here. 

w 3 j AS ° PPOSe , d t0 the COSmic ray eiectr °ns, which produce lower energy 7-rays by brehmsstrahlung. 
We do not explicitly consider the electron component here, but Prantzos and Casse argue that it should 
scale with the nuclear component as it does today, and using this scaling they calculate the contribution 
to the isotropic 7-ray background. 
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8 Conclusions 

In this paper we showed in detail that there are many uncertainties about early cosmic 
ray behavior and therefore a range of possible assumptions from which one must choose 
in building a model for this behavior. Furthermore, the variance in possible early cosmic 
ray models can have a significant effect of on the calculation of cosmic ray production of 
LiBeB (and 7-rays) in the early galaxy. The variance in the models derives both from 
uncertainties about present cosmic rays and from uncertainties about the evolution of 
their injection and acceleration, as well as confinement. 

Because these uncertainties are difficult to resolve, we conclude that future work in 
modeling early cosmic rays must allow for cosmic ray model variances and evolution. 
As shown m section 6.2, the choice of source spectral type is significant yet not settled 
even for today’s cosmic rays, with power laws both in total energy per nucleon and in 
momentum being allowed. Also, we show that cosmic ray evolution derives not only 
from changes in source strength, but also from changes in confinement, as detailed in 
sections 4 and 6.1. Evolution of confinement can have important effects, for example in 
altering the naively expected Be- to-Fe slope, as outlined here and shown in the detailed 
model of PCV. Additional considerations are noted in section 6.3, most notably the 
effects of inelastic nuclear collisions on confinement, as have been noted by Malanev and 
Butler (1993). 

In light of our findings we now turn to the work to date on early cosmic ray production 
of LiBeB, namely the relatively model-independent work of WSSOF and its sequels, and 
the detailed cosmic ray and chemical evolution model of PCV. The flatness of the curves 
in figure 13 justifies the claim in WSSOF that the B/Be ratio in cosmic ray models is 
indeed insensitive to the details of cosmic ray models. In addition, because Be and B are 
pure spallation (as opposed to a + a fusion) products, their production ratio dB/dBe is 
independent of [Fe/H] and varies only with the C:N:0 ratios of targets. To the extent 
that this ratio remains fairly constant, we can conclude that the B/Be ratio is insensitive 
to chemical evolution as well, and thus. We consequently find B/Be to be well-chosen 
as a model-independent signature of the spallation process. 

The Li isotopic ratio, discussed in SFOSW, 6 Li/ 7 Li shows a similar independence to 
cosmic ray model features but does have a moderate sensitivity to chemical evolution 
via the ISM ratio He/CNO, which sets the relative contributions of the a + a fusion 
versus spallation contributions to 6 Li/ 7 Li. We find that cTU/cPU is always within 
~ 30% of umt y for a11 spectra and escape parameters considered here. This uncertainty 
is much smaller than that of the 6 Li/ 7 Li ratio of Smith et al (1992). Consequently we 

find SFOSW’s assumption of a temporally constant d?Li/<FLi production ratio to be 
appropriate. 

However, we find that the “zeroth order model” of WS, SFOSW, and SFOSW to be 
overly simplistic not only in its lack of chemical evolution, but also in its lack of any 
kind of cosmic ray evolution. The effect on predictions of this model is unimportant for 
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the B/Be limits, but their calculations of the Li/Be ratio are, as we have noted, only a 
lower bound to the actual ratios. 

The PCV model does allow for chemical as well as cosmic ray evolution, and in 
particular they allow for evolution of cosmic ray confinement. We thus find their model 
indeed to be a useful examination of LiBeB production. Our caution is that their model 
is a specific one made with particular, albeit reasonable, assumptions. Their specific 
cosmic ray source model we find to be well-motivated and reasonable, but not a unique 
choice and we would urge that future work consider different momentum source indices, 
as well as source spectra in total energy. Furthermore, PCV admit that their particular 
implementation is to some degree arbitrary. In the first place, it is not obvious which 
parameters control A. Secondly, even given a particular scaling prescription, it is not 
trivial to model the input parameters (e.g. galactic scale height and magnetic field) 
accurately. PCV assume A depends upon the scale height of the collapsing disk, but 
they do not compute the collapse explicitly, and so impose an ad hoc time dependence. 
Clearly a firm model for A evolution is crucial, but at present is unavailable. We agree 
with PCV that more theoretical work in this area is needed. 

In summary, we feel that improved data on PopII LiBeB abundances with consistent 
controls on systematics might lead to an ability to constrain models for the origin and 
propagation of cosmic rays in the early galaxy. 
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FIGURE CAPTIONS 


Figure 1: (a) The cosmic ray proton flux and a spectra <f> compared to source spectra 
(<?(p) oc p ) scaled as A q. Note that in the limit of no energy losses (a good 
approximation at high energy) we have <j> = Aq , and so departure from this scaling 
indicates the effect of energy losses. Note also that the a flux scales very well with 
the proton flux as the He cosmic ray abundance 0.08. That this proportionality is 
accurate over the whole energy range reflects the similar energy loss behavior due 
to the scaling A/Z 2 . 

(b) As in ( a ), for q(T) oc (T + m p )“ 2 * 7 . 

Figure 2. (a) The transported cosmic ray proton flux for a range of confinement pa- 
rameters: A = 10, 30, 100, 300, and 1000 gem" 2 . The source is a power law in 
momentum, q(p) oc p~ 2 . 

(b) As in (a), for a source power law in total energy, q(T) oc (T -f- m p )“ 2 - 7 . 

Figure 3: The Li to Be production ratio as a function of pathlength. The solid line is 
for a source spectrum that is a power law in momentum, q(p) oc p~ 2 ; the dashed 
line is for a source power law in total energy, q(T) oc (T -fi m')~ 2 - 7 . We put fC /HI 
= [N/H] = -2.5, [O/H] = -2.0. 1 ' J 

Figure 4: The B to Be ratio as a function of pathlength, plotted as in fig. 3. 

Figure 5 : (a) The production ratio of Li to Be as a function of metallicity [Fe/H], for a 
source q(p) oc p~ 2 , and for a range of conement parameters: A = 10, 30, 100, 300, 
and 1000 gcm~ 2 . See discussion in the text regarding comparison to observations! 

{b) As in (a), for a source q(T) oc (T + m p )~ 2 - 7 . 

Figure 6: A plot of the pathlength distribution of equation (29), chosen in Garcia- 
Munoz et al. (1987) to provide the best fit to the observed B/C ratio. 

Figure 7: (a) The dLi/dBe ratio for A (T) = A 0 (solid line), and for A (T) given by 
equation (29), with Ao as indicated on the figure (broken line). The source spec- 
^ is q{p) oc p" 2 , and ISM abundances are [C/H] = [N/H] = -2.5 and [O/H] = 

(b) As in (a), for a source spectrum q(T) oc (T + m p )~ 2 - 7 . 

Figure 8: (a) The dB/dBe ratio for A(T) = A 0 (solid line), and for A(T) given by equa- 
tion (29), with A 0 as indicated on the figure (broken line). The source spectrum 
is q(p) oc p~ 2 and ISM abundances are [C/H] = [N/H] = -2.5 and [O/H] = -2.0. 

(b) As in (a), for a source spectrum q(T) oc (T -j- m p )~ 2 - 7 . 
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Figure 9: (a) The dLi/dBe ratio for A = A MC , the usual case (solid line), and for A (T) 
given by eq. (27), (broken line), where we have plotted against A esc . The source 
spectrum is q(p) oc p~ 2 and ISM abundances are [C/H] = [N/H] = -2.5 and [0/H] 
— - 2 . 0 . 

(b) As in (a), for a source spectrum q(T) oc (T -f m p )~ 2 - 7 . 

Figure 10: The ratios of Li to Be production rates, shown as a function of source 
spectral index 7 for a source spectrum in momentum and in total energy. Note that 
the yields from the momentum spectrum are very sensitive to the index adopted, 
and the ratio of rates can become unacceptably high in this case. In contrast, the 
yields from the total energy spectrum are relatively insensitive to spectral index, 
and are all at acceptably low levels. 

Figure 11: As in fig. 10 for the ratio of B to Be rates; the same trends are seen. Note 
however that the B/Be ratio is negligibly sensitive to the total energy index, and 
is only mildly sensitive to that for the momentum law. 

Figure 12: (a) Curves of dLi/dBe , as a function of A, for different spectral indices 7. 
The source is q(p) oc p~ 2 . 

(6) As in (a), for q(T) oc (T + m p )” 2 - 7 . 

Figure 13: As in fig. 12, for the dB/dBe ratio. 

Figure 14: Similar to figure 5. We plot dLi/dBe for each source type at A = 100 g cm” 2 . 
The data, as described in the text, has been corrected assuming the minimal 
allowed primordial 7 Li production, [ 7 Li]j3B = 2.0 and thus (linearly) subtracting 
this amount from the Li vales of table 1. See discussion in the text regarding 
caveats in comparing the theoretical curves, which are upper bounds, to the data. 

Figure 15: The dLi/dBe ratio for different determinations of the daughter nucleus 
kinematics. The solid curve is calculated assuming that the daughter nucleus is 
stationary in the center of momentum frame of the parents; this is the standard 
used throughout this paper. The dashed curve is calculated assuming that, in the 
rest frame of the light parent, the daughter kinetic energy per nucleon is equal 
to that of the heavy parent. The source spectrum is q(p) oc p~ 2 for curves (a 
and (6, and q(T) oc (T + m p )~ 2 - 7 for curves (c and (d. Additionally, we put 
[C/H]=[N/H]=-2.5, [O/H]=-2.0. 

Figure 16: (a) The dLi/dBe ratio for (solid line) cosmic ray CNO abundances equal 
to ISM abundances (pcno ~ Vcwoi the usual case), and for (dashed line) cosmic 
ray CNO abundances fixed to be solar (yg§ 0 > yg%). The source spectrum is 
q(p) oc p 2 and ISM abundances are [C/H] = [N/H] = -2.5, and [0/H] = -2.0 

(6) As in (a), for a source spectrum q(T) oc (T -f m p )~ 2 - 7 . 
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Figure 17: The dlA/ dBe ratio with incomplete LiBeB trapping, i.e. with 5'LiBeB = 
exp(j?/A) (solid line), and with complete LiBeB trapping (-SLiBeB = 1). The 
source spectra are as labeled, and we put and [C/H]=[N/H]=-2.5, [O/H]=-2.0. 

Figure 18: (a) The dhi/ dBe ratio in the absence (solid line) and presence (broken line) 
of a high energy plateau in the a + a — > 6,7 Li cross section. The source spectrum 
is q(p) oc p ~ 2 and ISM abundances are [C/HJ = [N/H] = -2.5 and [0/H] = -2.0. 
Note the small (~ 10%) increase with the plateau, 

(6) As in (a), for a source spectrum q(T) oc (T+m p )~ 2 - 7 . Here the difference 
is larger (a factor of ~ 2), as the flux is harder than that from the momentum 
flux, and so the high energy cross section behavior is more important. 
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